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Abstract 

The paper presents a model of a dynamic crack with a wavy surface. 
So far, theoretical analysis of crack front waves has been performed 
only for in-plane perturbations of the crack front. In the present pa- 
per, generalisation is given to a more general three-dimensional per- 
turbation, and equations that govern corrugation crack front waves are 
derived and analysed. 
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1 Introduction 

The paper analyses singular fields around a dynamic crack whose surface is 
slightly perturbed from the original plane configuration. Crack front waves 
in the plane of the crack were discovered numerically by Morrissey and 
Rice in [1], and later confirmed analytically by Ramanathan and Fisher [2], 
using the results of the perturbation analysis of Willis and Movchan [3]. 
Experimental observations of persistent crack front waves were reported by 
Sharon, Cohen and Fineberg [4]. The more general development of Willis 
and Movchan [5] and Woolfries et al. [6] extended the analysis to a crack 
propagating through a viscoelastic medium. The perturbation formulae for 
the stress intensity factors, specialised to a plane strain formulation, have 
been used by Obrezanova et al. [7] in the stability analysis of rectilinear 
propagation. A quasi-static advance of a tunnel crack under a mixed mode 
loading has been analysed by Lazarus and Leblond [8]. 

The aim of the present paper is to develop a model describing corrugation 
(out-of-plane) waves along the front of a moving crack. This work is based 
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on the ideas of the earlier publication by Willis [9] . The plan of the paper 
is as follows. We begin, in Section 2, with the description of the geometry, 
governing equations and perturbation functions. A summary of the first- 
order approximations for the stress intensity factors is presented in Section 
2.2. Section 3 includes the study of the corrugation waves in the first-order 
asymptotic approximation for a basic Mode I loading. In Section 4, we 
derive the dispersion equation for crack front waves in the mixed mode I- 
III loading. The technical appendix contains an outline of the fundamental 
integral identity, and the expressions for effective tractions. 

2 Basic perturbation formulae 

For a linearly elastic medium, we consider a semi- infinite crack with a slightly 
perturbed surface. The unperturbed configuration of the crack at time t is 
defined by 

So(t) = {x : — oo < x\ < Vt, — oo < %2 < oo, X3 = 0}, (1) 

where V is a constant crack speed, which does not exceed the Rayleigh wave 
speed. The perturbation is introduced through deviations of the crack front 
in both in-plane and out-of-plane directions. The perturbed surface of the 
crack at time t is 

S £ (t) = {x : —00 < x\ < Vt + £(p(x2,i), 

— OO < X2 < OO, X3 = £1p(xi — Vt,X2)}- (2) 

The functions tp and ip are smooth and bounded, and e is a small non- 
dimensional parameter, < e < 1. It is helpful to use the moving- frame 
coordinates, so that X = x\ — Vt. 

It is assumed that the medium is loaded so that a stress er nc and a 
displacement u nc would be generated in the absence of the crack. The crack 
induces additional fields cr, u. They satisfy the equations of motion and the 
traction boundary conditions on the crack faces: 

a ij,j ~ P^i = 0) i = 1) 2, 3, outside the crack (3) 

and 

TijUj + afjUj = 0, on the crack faces, (4) 
and correspond to waves outgoing from the crack as X3 — > ±00. 
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2.1 Local coordinates and asymptotics for stresses 

At a point x° = (x^x^x®), which is on the crack edge at time t, so that 

x\ = Vt + e<p{x° 2 ,t), xl = e${x\-Vt,xl), 
we define a coordinate system such that 

3 

x — x y ] x i e ii 
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where 



Here ^* denotes ip evaluated for x\ = Vt. The above transformation involves 

a shift to the crack edge and a further rotation of coordinate axes. 

In the new frame, the stress components (u' i3 ) have the asymptotic form 

a' i3 (x[, 4, 0) ~ (K^ + sK^)/(2nx[) 1 / 2 - (p/ 0) + AP,) 

+ (Af ) + A^(x;) 1 / 2 -(^ (0) +AF,y i + (iv( ) + AiV i )(x , 1 ) 3 / 2 , i = 1,2,3. (6) 

The first-order asymptotic approximation of stress-intensity factors was 
constructed and studied in [3], [10], [11], [12]. In Appendix we include a 
description of the fundamental identity, which is essential for this work. We 
also require the dynamic crack face weight function [U], as defined in Ap- 
pendix. The field [U] has a singularity proportional to X~ 1 / 2 H(X)5(x2)5(t) 
as X -> 0. 

2.2 First-order perturbations of the stress intensity factors 

We begin with the first-order approximation for the stress intensity factors, 
when 

Kj ~ Kf + eKf ] ,j = I, II, III. (7) 

For the Mode-I unperturbed case, Kjf = Kj°jj = 0, and the perturbation 
terms are defined by (see [11], [12]) 
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+ [t^lii * + Phi * (P 2 (1) ) - (C/)3i * [P^L (8) 

+ [U] 12 * (P^) + [U} 22 * (P 2 (1) ) - ([7} 32 * [P^], (9) 

-(C/}23*[P 2 (1) ] + ^]33*(P 3 (1) )- (10) 

The matrix Q is a block-diagonal matrix defined in [11]; other functions 
that appear in the above equations are 

013 = En + ^Sl2, wis = T^t 1 + + 2 £) " 2 ' 

-23 = ^(l + /3 2 )(« 2 -/3 2 )-l, 
v 4a/?-(l + 2a 2 -/3 2 )(l + /3 2 ) -2(1 + fi 2 - 2a/?) 

a 2 = 1 - V 2 /a 2 , /3 2 = 1 - y 2 /6 2 , i?(F) = 4a/? - (1 + /3 2 ) 2 . (11) 

Here, a and 6 denote the speeds of longitudinal and shear waves, respectively. 
The representations for the effective tractions P^\i = 1,2,3, are given in 
Appendix. 



2.3 Crack front waves confined to the plane x 3 = 

Assuming that the out-of-plane deflection is not present (tp = 0) , we consider 
a first-order in-plane perturbation of the crack front and loading in Mode I, 
so that 0"i3 = <72 3 = on the plane x 3 = 0. In this special case, the only 
non-zero stress intensity factor is Kj, and the corresponding perturbation 
formula reduces to 

According to the Griffith energy balance equation, the energy flux Q into 
the crack edge is constant, denoted here by Q c : 

G = l -^h{v)K 2 = g c . (13) 
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Here, v is the local crack speed (to the first-order approximation, v = V+eip) 
and fi(v) is a known function (e.g., [13]): 

ip' a(v) 

A <"> W (14 » 

Expanding the Griffith energy balance equation (13) to order e, we obtain 

f'(V) 

2Q 33 * V + 7^7^ + 2rmp = 0, (15) 

where m = (it /2) 1 / 2 A 3 °^ j Kj°\ Applying the Fourier transform with re- 
spect to t and X2 we deduce that a non-zero solution is possible only if the 
dispersion relation 

f'(V) 

2Q 33 (oj, k) - M-f-!- + 2m = (16) 

Ji\Y) 

is satisfied. Here, the Fourier transform Q 33 is a homogeneous function of 
degree 1 in (uj, k). At high frequency and large wavenumber, the third term 
in the above equation can be neglected. Such an equation can be solved for 
oj/k, and a real root represents a speed of wave propagating along the crack 
front. This computation was performed by Ramanathan and Fisher [2]. 



3 Corrugation waves for a Mode-I basic loading. 
First-order analysis. 

Can a Mode-I basic loading generate a corrugation wave propagating along 
the crack front? This case corresponds to a non-zero out-of-plane pertur- 
bation characterised by the function ip(xi — Vt,X2)- Crack stability with 
respect to out-of-plane deflections can be studied, once a fracture criterion 
is identified. 

If we suppose that Kjj = then, to lowest order, ip must satisfy 
Kflty) = 0, where K$ is given by (8). The proposition that the crack 
propagates so as to maintain Kjj = together with the Griffith energy bal- 
ance has recently received theoretical support, on the basis of a version of 
Hamilton's principle [14]. 

Assuming that the in-plane perturbation of the crack front equals zero, 
we look into stability against out-of-plane deflections. It is also assumed 
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Figure 1: Level curves of the function W = |Qn6i3 — i(a;/V)a;i3|, for V/b = 
0.69 and v = 0.3. 

that uj = k\V and &2 are large. The leading-order approximation of the 
stress intensity factor Kn yields 

K$ = {-QnQia + i(w/V>i3}tf}V = 0. (17) 

This relation is homogeneous of degree 1 in w and &2, and so is non- 
dispersive. 

The numerical study of equation (17) produced the following results. 

• For crack speeds V greater than a critical value V c (which is close to 0.6 
of the Rayleigh wave speed) there is a value i] = u/\k2\ with small, negative, 
imaginary part that satisfies (17). The position of the root is shown in 
Figure 1; the calculation is produced for the case of V/b = 0.69, and the 
diagram shows the level curves of the modulus of the expression in the curly 
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brackets on the left side of (17). 

Figure 1 is accompanied by a three dimensional surface plot, shown in 
Figure 2, of the function W = \Qu@i3 — i(co /V)cois\; the surface touches 
the ry-plane at the point corresponding to the root of equation (17). 




Figure 2: Surface plot of the function W = |Qn@i3 — i(w/V)wi3|, for V/b = 
0.69 and v = 0.3. 

• The "corrugation wave" suffers slow attenuation as it propagates. The 
imaginary part of 77, which characterises the rate of attenuation of the "cor- 
rugation wave", is shown in Fig. 3 for different values of the crack front 
velocity V, and it decreases with V. 
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Figure 3: The imaginary part of rj as a function of the crack front velocity 
V. 

4 First-order coupling between in-plane and out- 
of-plane crack front perturbations for mixed Mode 
I-III loading 

Here, we assume that Kjf = 0, whereas and Kj°jj are non-zero for a 
half-plane crack propagating with constant speed V (unperturbed configu- 
ration). To first order, the stress intensity factors are represented by the 
formulae (7), where the perturbation terms j = I, II, III, are defined 

by (see [11], [12]) 

K$ = -Qn * r&i3Kf - V>i3*f } - r [sn + ^n) 4° } /| 
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+[U] U * (P^) + [U] 21 * (P 2 {1) ) - (U) 31 * [P 3 (1) ], (18) 
^ffl = -Q12 * r^isK^ - ^> 23 *f ) + Q22 * (<pKffj) + 

+ [C/]i2 * (P^) + [^22 * (P 2 (1) ) - (U) 32 * [P 3 (1) ], (19) 

^ = Q33 * ^} 0) + (|) 1/2 ^4 0) - ^(1 - S Si2 K } \/I 

-2^ 2 tf$ - (£/} 13 * [P x (1) ] - (C/} 23 * [P 2 (1) ] + [U} 33 * (P 3 {1) >. (20) 

We shall use the criterion of local symmetry Ku = 0, together with the 
Griffith energy balance equation 

G = (2 f i)- 1 f I (v)K 2 + (2fir 1 f III (v)K] II = Q C = const. (21) 

Taking into account that, to first order, v ~ V + £92, we deduce 

= (2^)- 1 //(V^)(< ) ) 2 + (2^)- 1 /m(^)(^° / ) / ) 2 

+e(2^- 1 ^f' I (V)(KPf + 2/ / (F)k} 0) k} 1) 

+^ jj (F)(k}? j ) 2 + 2f III (y)Kj t %K& I ) + 0(e 2 ). (22) 
It follows from (21), (22) and the local symmetry criterion Kn = that 

+2f in {V)Kf} I K$ I {ip, V) = 0, (23) 

tf£W)=0. (24) 

The above equations define the coupling between the in-plane and out-of- 
plane perturbations of the crack front. 

Applying the Fourier transform with respect to t and x 2 and assuming 
that oj = k\V and k 2 are large, we deduce 

Tp{ (2f I (V)Q 33 - kjfKvj) (KP) 2 + (2f IH (V)Q 22 - kjf' in (vj) (K^) 2 } 
+rKf ) K { $ I {2f III {V) (-Q 12 0i 3 + ik 2 oj 23 ) + 4ik 2 f T (V)} = 0, (25) 
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{-Qn&i3 + i(u/V)LJ 13 }ip*K\ Uj + (Q 21 + ifc)?^ = 0. (26) 

The system (25), (26) is linear in Tp and ip*, and it possesses a nontrivial 
solution if and only if the matrix of this system is degenerate. This yields 
the following dispersion relation: 

{-Qnei3+i(w/y)^i3}{2//(F)Q 3 3-ia;/KV r )+(2////(y)Q 22 - iuf' ni (V)) IC 2 } 

-JC 2 (Q 21 + ik 2 ) {2f HI (V) (-Q 12 e 13 + ik 2 u 23 ) + 4ik 2 fj(V)} = 0. (27) 

Here /Co = Kfii/K-f^ ■ The above dispersion equation, connecting co and k 2 , 
is to be analysed numerically to identify possible crack front waves associated 
with the external mixed mode I-III load. 

Appendix. Fundamental identity and effective trac- 
tions. 

Here, we briefly describe the method developed in [3], [10], [11]. We use the 
relation 

u = -G*<r, (Al) 

where u and <x denote the values of the displacement vector (m) and the 
traction vector (0^3) on the surface x% = of the half-space X3 > 0; G is the 
Green's matrix function. The symbol * denotes convolution over x±, x 2 and 
t. It is assumed that all waves emanate from the surface X3 = 0. A similar 
identity applies to the half-space X3 < 0, with G being replaced by — G T . 

Three column vectors like u can be written side by side to form a matrix 
U(+0), and similarly Xl(-|-0) represents the matrix formed from the three 
corresponding vectors a. Then 

U(+0) = -G*S(+0). (A2) 

The argument (+0) signifies values on the boundary of the upper half-space. 
Applying similar reasoning to the identity for the lower half-space X3 < 
gives 

U(-0) = G T *S(-0). (A3) 

Next, we note that 

{U(+0)} T * <t(-0) = -{£(+0)} T * G T * <t(-0) 
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= -{E(+0)} J *u(-0), (A4) 
{U(-0)} T * <t(+0) = {S(-0)} T * G T * <t(+0) 

= -{£(-0)} T *u(+0). (A5) 

Subtracting the second line from the first and rearranging gives the identity 

[U] T * (a) - (U) T * [a] = - [S] T * (u) + (S) T * [u] , (A6) 

where (/} = i(/(+0) + /(-0)) and [/] = /(+0) - /(-0). 

In the moving frame associated with the crack edge, we use the coordi- 
nate X = x\ — Vt. The operation of convolution survives, with functions 
regarded as functions of X, x-i , t and the convolutions taken over these new 
variables. 

For the unperturbed crack problem, 

[a] = 0, [u] = when X > 0, a = (a) = -er Ilc when X < 0. (A7) 

We interpret equation (A6) relative to the moving frame, and perform fac- 
torizations of the Green's function so that U and 5] display the related 
properties 

[S] = 0, [U] = when X < 0, E = (£) = when X > 0. (A8) 

Equations (A2), (A3) yield 

[U] = -(G + G T ) * (S), (U) = -§(G - G T ) * (S). (A9) 

The first of these relations defines a Wiener-Hopf problem; the second then 
gives (U) directly. The Wiener-Hopf problem uncouples into two sub- 
problems. One, associated with the opening mode I of the crack, is a scalar 
problem. It was solved in the case of elasticity in [3], and for a viscoelastic 
medium in [6]. The remaining problem involves modes II and III, coupled. 
It was solved in [10]. 

The field [U] has a singularity proportional to X~ 1 / 2 H(X)5(x2)S(t) as 
X —> 0. With the constant of proportionality chosen as (2/ir) l / 2 l, we call 
[U] the dynamic weight function for the crack problem. With this choice, 
letting X — > +0 in the identity (A6) generates 

K = x lim {<Uf * [<x(°)] - [Vf * <<x(°)>}, (A10) 
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where K denotes the vector of stress-intensity factors (Kjj, Kjjj, Kj) T . 
The matrix function (U) represents a dynamical version of Bueckner's non- 
symmetric weight function, as described in [15] and [10]. 

We assume that the unperturbed steady-state crack is subjected to a 
Mode-I loading, and the unperturbed displacement field is a vector function 
u (o) _ u (°)(xi — Vt,X2,xs). We can write the resulting displacement field 
in the form 

u-uW+euW 

where e is a perturbation parameter. 

The effective tractions P^ := —(Tja{vft))\x z =a, i = 1,2, 3, have the form 
(see formula (4.11) of [11]) 

k=l 
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